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1 Introduction 

Imagine three-space to be filled with an elastic material and let rj(t) : R 3 — > R 3 
be the map which takes each material point from its position at time zero 
to its position at time t. Then the curve of maps, t i— > rj(t), will completely 
describe the motion of the material. 

We shall assume that the material is incompressible, which in mathemat- 
ical terms means that r](t) must satisfy the equation: 



where J means the Jacobian determinant. We shall further assume that the 
material is "neo-Hookean" or that when deformed by a map rj : R 3 — > R 3 , 
it acquires potential energy: 



V(v) = ~ {dirfdirf — 3) (repeated indices summed from 1 to 3) (2) 



J(v(t)) = 1 



(1) 




From ([[J) - (|2]) one can derive the equations of motion which rj(t) must 
satisfy, (see @). They are: 



V\t)(x) = (A V \t))(x) + (d iP (t))(r](t)(x)) (3) 
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where " • " means time derivative, "A" means Laplacian, and p(t) is a function 
on R 3 which is chosen so that (1) will hold, (as we shall see below). 

It is known || that given initial data 77(0) and r](0), there is a positive 
T and a unique curve r](t) (0 < t < T) satisfying (1) and (3).[] In this 
announcement, we shall show that if 77(0) is sufficiently near the identity 
map and if fj(0) is small, then one can take T = 00. Thus if the initial 
deformations are small, the equations have a unique solution for all time. 

To provide context to this theorem we now discuss some related results, 
each of whose proofs contains ideas that we will use for our proof below. 
First we consider the case of compressible elastic media. If the potential 
energy is as in (|j), the resulting equation is the (linear) wave equation, so 
of course solutions persist for all time. However if V{rj) is not quadratic, the 
resulting equation is not linear, so it is possible that solutions will develop 
singularities. F. John |J M analyzes this situation in considerable detail for 
homogeneous isotropic materials. He finds that for this case the situation 
is essentially the same as for the simpler case of a single quasi-linear wave 
equation, which has been studied over many years by Klainerman, John, 
et al. (see [0] for a survey). Thus we proceed with a brief discussion of 
quasi-linear wave equations (some of which we shall need for the proof of our 
theorem) . 

Let / = f(x°, x 1 , x 2 , x 3 ) (identifying t with x°) and let /' and /" denote 
respectively its first and second order partial derivatives. Let / satisfy: 

□/ = F(f, f") (4) 
where □ = 8% — 8f — 8\ — <9f , F is linear in /", and: 

F{f',f") = 0{\f'\{\f'\ + \f"\)) 

near zero. Then John and Klainerman prove the following: || 

Theorem 1. Let f°(x) and J 1 (x) be smooth functions on R 3 with 
compact support. Then there exist positive constants A, B, and Eq such 
that if e < Sq, then equation (|j) with initial data f(0,x) = ef°(x) and 
dtf(0, x) = ef'(x) has a unique solution for < t < T with T > Ae B ^ £ . 

Thus a quasi-linear wave equation has solutions for a very long time if 
the data are small. Moreover as John points out, (0 p. 27) this result is 

"solution" we will always mean smooth or classical solution. We do not consider 
weak solutions in this paper. 
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sharp. That is letting: 

and choosing f 1 with support in the unit ball, we find that (§]) does not have 
a smooth solution on [0, T) for T > e l ^ L where: 

Theorem 1 is proven by looking at various weighted spatial L 2 norms 
of / and its derivatives and showing that they obey certain differential in- 
equalities. In particular one shows that if N(t) 2 is a certain sum of terms of 
type: 

/ \(x a )\dp) k f{t,x)\ 2 dx a,P = 0,1,2,3, t < k (5) 

with t identified with x°, and ), then one gets an inequality of 

type: 

^-N{t) < const. (1 + t)- l N(t) 2 (6) 

Using this and the fact that ^(l + t)^ 1 ^ = log(l + T) one gets the necessary 
estimates to prove Theorem 1. 

On the other hand, the sharpness of Theorem 1, or the necessary devel- 
opment of singularities is shown by looking at spherical means of functions 
depending on / and its derivatives. Using these one constructs a function 
W{t) which obeys: 

j t W{t) > const. (1 + t)- l W{tf (the reverse of (6)) (7) 

and from this gets an upper bound on the time of existence of a smooth 
solution. 

John shows that for compressible materials, the situation is the same as 
for solutions of (|4|). He uses the displacement u(t, x) = r](t)(x) — x and notes 
that it satisfies an equation: 

d 2 u = dAu + c 2 Vdiv(w) + F(u', u") (8) 
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where "V" means gradient, "div" means divergence, F is like the correspond- 
ing function in (£|) and C\ and C2 are positive constants. For (||) one can get 
estimates like those of (Q) and thus find that given initial data of compact 
support, solutions exist as in Theorem 1. (see 0]) Here too John shows that 
the result is sharp, (see |J) Taking spherically symmetric data one gets a 
solution of the form: 

u(t, x) = (p(t, \x\)x (9) 

where ip(t, \x\) satisfies a quasi-linear wave equation. Certain combinations of 
second derivatives of ip obey an inequality like (^) and thus ip{t, \x\) remains 
smooth for only a finite time. 

For incompressible elastic materials all results that we know of are for a 
finite time which depends on the data. proves such a result for materials 
filling R n , and [El and || show a similar result for initial-boundary value 
problems. Furthermore, we do have examples where singularities develop. 
(||, section 7) 

Another relevant case is that of incompressible fluids (the case V{rj) = 0) 
for which there is a much larger literature. First there is the old result of 
Wolibner fllU| . It uses estimates based on conservation of vorticity to get 
existence of solutions for all time for fluids in R 2 . The same principle is used 
by Kato |5| and Yudovic [ITJ who prove existence for all time for fluids in 



two-dimensional domains with boundary. This analysis is further extended 
by Beale, Kato and Majda |TJ who show that existence of solutions for all 
time is guaranteed in R 3 provided the vorticity of the flow obeys a certain 
bound. 

The proof of our result uses a combination of all of the above techniques. 
We first consider the equation in terms of the displacement. Applying the 
curl (or exterior derivative) to that, we get a quantity somewhat analogous 
to the vorticity of a fluid, which satisfies an equation similar to @. Using 
norms like N(t) above we get an inequality stronger than and from that 
we get existence of solutions for all time. 

We would like to generalize our proof so that it covers incompressible 
materials which are not neo-Hookean, but have not been able to do so as yet. 
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2 Definitions and Statement of the Theorem 



We shall be concerned with functions and vector fields on R 3 , and we shall 
use the exterior derivative, "cT , on either of these. Thus for a function p, dp 
is the same as Vp and for a vector field u, du is the same as the curl of u. 
Also we will denote by 8 the formal adjoint of d so that 8u is the negative of 
the divergence of u. 

Furthermore if L is any differential operator and t] : R 3 — > R 3 is a 
diffeomorphism, we shall define L v by: 

L v f = (L(for ] - 1 ))or 1 . (10) 

// r)(t) is a curve of diffeomorphisms, a direct calculation yields: 

j t (L v{t) f) = [vV,L] v{t) f (11) 

where [ , ] denotes the commutator and v is the vector-field defined by 

fi(t)(x)=v(ri(t)(x)). (12) 

Finally for functions on space-time we define the operators: 

Ti = di i = • • • 3 
T 3+ j = x k di — x l dk where {ijk} is an even permutation of {123} 
r 6+fe = td k + x k d fc = l---3 

3 

Tio = td Q + J2 xi 9i 
i=i 

Also for any multi-index a = (a ■ ■ ■ a w ) we define the higher order op- 
erator: 

T a = r^T" 1 •••r° 10 . (13) 

Tq to r 3 of course commute with the wave operator □, and since T4 to 
Tg are infinitesimal generators of the Lorentz group, they commute with □ 
also. However Tio generates similarity transformations and for it we have 
the identity: 

[r ,n] = -2n. (14) 
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Using the T p we define certain 1? norms on functions of space time which 
depend on time and are stronger than the Sobolev norms. We write: 



\Us= E [ \r a f(t, x )\ 2 dx (is) 



\a\<s 

(LS stands for Lorentz-similarity) . We shall also need supremum norms for 
functions and their derivatives and for this we write: 

\\f\\ ck = £ sup \D a f(x)\ 

\a\<k l6R3 

where D a = d^d^d^. 

Having introduced the necessary notation, we can now discuss our ba- 
sic equation @. We introduce q(t)(x) = p(t)(rj(t)(x)) and then write the 
equation as: 

fj = An + d v q. (16) 

We would like to show that q is determined by rj and its derivatives, and to 
do so we first take a time derivative of equation (1). This gives: 

$v(v) = 0. (17) 
Then applying a second time derivative and using flTT| ) gives: 

^(77)-[v.V,<y]^ = 0. (18) 
where v is defined as in (|r2|). Also a direct calculation yields: 

[v ■ V, 5] v i) = tr((D r? ?)) 2 ) where "tr" means trace of a matrix. (19) 
Now applying 5 V to (]T^) and using (|TJ) and flT9"|) we find: 

tr(( J D r ,77) 2 ) = ^A ?7 -A ?) g. (20) 



Thus given r)(t) and r)(t) we get q(t) as the solution of (p0|) and it follows 
that the right side of (|nj) is a function of r\ and rj. Our main theorem is: 

Theorem 2: Given r/ and T]i smooth vector fields on R 3 such that r] — Id 
and rji are rapidly decreasing; there exists a maximal T > and a unique 
curve rj(t) (0 < t < T) which satisfies ([[]) and ( fT6|) (g defined as in (pOD) such 
that 77(0) = r]o and ?)(0) = r/i. Furthermore there exists e > such that if 
\\t]q — Id\\ LS ,5 and ||??i||ls,4 are less than e, then T = 00. 
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3 Outline of Proof 



We let u(t, x) = rj(t)(x) — x be the displacement as denned in §1. Then (|16|) 
is equivalent to: 

□u = d^q. (21) 

Also since t](t)(x) = x + u(t,x), we find that the equation (1) is equivalent 
to: 

8(u) = Q(Du) + det(Du) (22) 

where 

3 / du 1 d u i+1 \ 
Q(Du) = y^det[ n * , l , (indices taken mod 3) (23) 

fri V di + lU w u J 

and "det" means determinant. Thus (^1]),(^) is equivalent to ([!]), (|3]) and 
from we have the following short time existence theorem. 

Theorem 3: Given u° G # 4 (R 3 ,R 3 ) and u 1 G if 3 (R 3 ,R 3 ) such that 
5(u°) = Q(Du°) + J(u°) and ^(tt 1 ) = ; there exists T > and a unique solu- 
tion of u(t, x) forO<t<T such that u(0, x) = u°(x) d t u(0, x) = 
v}{x). Furthermore if u° G if s+1 (R 3 ,R 3 ) and u 1 G iJ s (R 3 ,R 3 ) then u G 
n| =0 C fc ([0,T), ^ +1 - fe (R 3 ,R 3 )). 

For our estimates we will need the || \\ls norms, which are examples of the 
norms N(t) of §1. Hence we need to know that u(t,x) goes to zero rapidly 
as \x\ — > oo. Thus we introduce H s,e norms: 

t = f \D s f(x)\ 2 (l + \x\ 2 ) e dx. (24) 

A slight modification of the proof of [0] shows that in Theorem 3, H s can be 
replaced by H s,i for any positive integer I. Thus norms of the form N(t) will 
be finite for solutions u(t,x). 

To prove the last statement of Theorem 2 it suffices to show that u remains 
bounded as t gets large. To do this we will derive an equivalent equation for 
du and show that solutions of that equation must remain bounded. 

To derive an equation for du we first apply d, q to (^T|) getting: 

d v a u = 0. (25) 
However, by direct calculation we find that for any v : R 3 — > R 3 , 

d v v = A((Drf)-^Dv) (26) 
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where (Drj) 1 denotes the matrix inverse of Drj and A applied to a matrix 
means the matrix minus its transpose. Thus fl25|) can be written: 

A((Dri)-*UDu) = 0. (27) 

Since Dr/ — Id + Du (and therefore (Drf)' 1 — Id = —Du(Drj)^ 1 ) we can 
rewrite (^7j) as: 

A{UDu) = A{Du{Dr])- l UDu) 

or 

Udu = A(Du(Dr])- l nDu) (28) 

Using LS norms we will derive estimates for (p8|) which show that it has 
solutions for all time. 

To analyze ( p8|) we would like to eliminate the highest order terms on the 
right side - namely UDu, and to do this must first express Du as a function 
of du. Letting uj = du and using (B2), we find that u satisfies the elliptic 
system: 

du = u, 5u = Q{Du) + det(Du). (29) 

Proposition 4. There exists a ball B about zero in i7 3 (R 3 ,R 3 ) such 
that given u G i^ 2 (R 3 ,R 3 ) there is at most one u in B which satisfies $2$). 

Proof. Use the fact that u is determined by du and 5u plus the fact that 
Q(Du) and det(Du) are 0(\Du\ 2 ). 

Proposition 4 defines a map G{oS) = Du where u is the solution of (|2H|). 
Thus (^) can be written: 

Udu = A{Du{Dri)- l UG{du)). (30) 

We proceed to analyze OG(du), our first step being to compute: 

d a G(du), a = 0, ... ,3 

Applying d a to (^) we get: 

dd a u = d a u 5d a u = Q'(Du)Dd a u + det'(Du)Dd a u (31) 

where Q' and det' are the derivatives of Q and det. Thus d Q G(u>) = Dd Q u 
where d a u is the solution of fl3~T|). We write Dd a u = G'{uo)d a u. 
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Applying d a again we find that d a u satisfies: 
dd 2 a u = d 2 a u 

5d 2 a u = Q\Du)Dd 2 a u + det'(Du)Dd 2 a u (32) 
+Q"(Du)(Dd a u, Dd a u) + det" \Du)(Dd a u, Dd a u) 

where Q" and det" are second derivatives of Q and det. Thus we can write 
d 2 u = c a + b a where c a solves: 

dc a = d 2 a u Sc a = Q'(Du)Dc a + det\Du)Dc a (33) 

and b a solves: 

db a = 

5b a = Q'(Du)Db a + det 1 \Du)Db a (34) 
+Q"(Du)(Dd a u, Dd a u) + det" (Du)(Dd a u, Dd a u). 

In particular Dc a = G'(u)d^u. Then summing over a we get: 

3 3 

UG{u) = Dico-^cJ + Dibo-^k) (35) 

i=l i=l 

= f DC + DB 

Here 

DC = G'(u)(ndu) (36) 
Also combining with ( p0|) we get: 

□du - A(Dn(Dr7)- 1 J DC) = A(Du(Dr})~ 1 DB) (37) 

so using ( |3~6"D we get: 

Udu - A{Du{Drj)- 1 G'{w)Udu) = A{Du{Dr 1 )- 1 DB). (38) 

The second term on the left side of (|38|) can be thought of as a linear 
operator applied to Ddu and if Du is small it will have small norm. Hence 
the map *f>(ndu) = Ddu — A(Du(Drf)~ l G' '(du)Odu) will be a linear map 
which is near the identity. Thus $ will be invertible and we can write: 

ndu = $- 1 (A(Du(Dr))- 1 DB). (39) 
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Our theorem can now be proven by getting estimates for the solution of 
©• 

We note that DB, which is a sum of Db a , depends on Dd a u in a bi-linear 
fashion. Therefore, letting dDu denote the four partials of Du, we get: 

Du{Drf)- l DB = 0(\Du\\dDu\ 2 ) (40) 

and since $ goes to the identity as Du — * 0, we find: 

^- 1 {A{Du{Drf)' 1 DB) = 0{\Du\\dDu\ 2 ) (41) 

also. 

Let E denote the energy associated with the wave equation, i.e. 

E{g) = \j^{\Dg\ 2 +\d t g\ 2 ). (42) 
Then f t E(g) = f K3 (n g )(d t g). Therefore: 



d 
di 



E(du)= J(^- 1 (A(Du(Dr])- 1 DB)),d t du). (43) 



Then applying T a to fl3T3| ) (with |a| < 4) and using ( |2 ) and the commutation 
properties of □ and T a , one gets: 



d 
dt 



E(T a du) = j \v a ^{A{Du{Dri)- 1 DB)),T a d t du) (44) 
< K\\d t du\\ LSA \\ddu\\ LSA ||Z>u||ca 



for some constant K. 

But Klainerman's inequality (see |7|], appendix 2) combined with the es- 
timates for Du = G(du) tells us that: 

\\Du\\ 2 c s < — ^— J2 E{T«du). (45) 

\a\<4 

Let E 4 {du) = J2 H < 4 E{T a du). Then E 4 (du) bounds \\ddu\\ 2 LS4 , so from (g) 
and (|45| ) we find that there is a constant K such that: 

^-E 4 (du) < K(l + ty 2 E 4 (du) 2 . (46) 
10 



But integrating this inequality we find that if E&(du)(0) < j^, then: 

1 - K£; 4 (dM)(0)(l - 

for all positive t. 

From this it follows that if du (and its derivatives) is small at time zero, 
then it must remain bounded for all time. Thus the solution of (|39|) remains 
bounded, and hence the solution of fl2"T| ) must remain bounded as well. 
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